Spin squeezing, entanglement and quantum metrology with Bose-Einstein 

condensates. 



o 






C0 



T3 

a 
o 
o 



> 

ON 

in 

CO 

in 

en 

o 



? 



Christian Gross 
Kirchhoff-Institut fur Physik, Universitdt 
Im Neuenheimer Feld 227, 69120 Heidelberg, German^ 

Squeezed states, a special kind of entangled states, are known as a useful resource for quantum 
metrology. In interferometric sensors they allow to overcome the "classical" projection noise limit 
stemming from the independent nature of the individual photons or atoms within the interferometer. 
Motivated by the potential impact on metrology as wells as by fundamental questions in the context 
of entanglement, a lot of theoretical and experimental effort has been made to study squeezed states. 
The first squeezed states useful for quantum enhanced metrology have been proposed and generated 
in quantum optics, where the squeezed variables are the coherences of the light field. In this tutorial 
we focus on spin squeezing in atomic systems. We give an introduction to its concepts and discuss 
its generation in Bose-Einstein condensates. We discuss in detail the experimental requirements 
necessary for the generation and direct detection of coherent spin squeezing. Two exemplary exper- 
iments demonstrating adiabatically prepared spin squeezing based on motional degrees of freedom 
and diabatically realized spin squeezing based on internal hyperfine degrees of freedom are discussed. 



I. INTRODUCTION. 

Interferometers, which are among the most precise 
measurement devices known today, are based on the in- 
terference of two coherent atomic or photonic modes. 
The quantity to be measured is mapped to the relative 
phase between the two modes, explicitly using the wave- 
like nature of the resources (atoms or photons). However, 
in the actual detection process their particle-like charac- 
ter becomes important; the population of the two modes 
can be at best be counted one by one. Traditional in- 
terferometers use uncorrelated resources in two coherent 
modes which results in a Poissonian distribution of the 
population difference at the interferometer output. The 
standard deviation of this distribution yN is the so called 
projection noise or shot noise limit. 

The best theoretically possible phase precision for un- 
correlated resources and linear interferometers is given 
by the classical Cramer-Rao bound of parameter estima- 
tion [TJ [2]. In the context of (classical) interferometric 
phase estimation this is known as the standard quan- 
tum limit. It can be reached with coherent states, the 
quantum states that are naturally generated by coherent 
splitting of a quantum mode, as it is routinely done by 
beam splitters in optics (with one input port unused). 
Indeed, modern interferometers measuring for example 
time [3] , position [4j [5] or magnetic fields [6] are limited 
by this quantum noise leading to a maximal interfero- 
metric phase precision of Acp = 1/y/N. 

Several strategies have been proposed to overcome this 
standard quantum limit for metrology by using quantum 
correlated, that is entangled, resources within the inter- 
ferometer [7] . One of these strategies is to use squeezed 



Present address: Max-Planck-Institut fur Quantenoptik, Hans- 
Kopfermann-Strasse 1, 85748 Garching, Germany; chris- 
tian. gross@mpg.mpq.de 



quantum states [8] which have been successfully gener- 
ated in quantum optics [9 . The analogue for atom inter- 
ferometry are spin squeezed states [TOl[TT]. Spin squeezed 
states can be obtained from coherent spin states by a 
continuous transformation which gradually redistributes 
the quantum noise among two orthogonal spin directions. 
This property makes them promising candidates for their 
experimental realization since the degree of squeezing can 
be balanced with the increased sensitivity to decoherence. 
Interferometric precision is increased to Acp = £s/V~N for 
spin squeezed states with a squeezing parameter £5. 



Many-body quantum mechanics sets an ultimate limit 
to the phase precision for general interferometric mea- 
surements. This is the so-called Heisenberg limit, for 
which an 1/N scaling of the interferometric precision can 
be achieved [7J. The Heisenberg limit can be reached 
with spin squeezed states in the limit of infinite squeez- 
ing where spin squeezed states become identical to twin 
Fock states. The potential gain in phase precision is 
enormous: In an atomic clock, for example, one mea- 
sures an energy difference Huj for a certain time t (2ith is 
Planck's constant). The "classical" uncertainty scales as 
Auo 2 = 1/tN, while Heisenberg limited metrology would 
allow for an error of Au 2 = 1/tN 2 . Assuming a fic- 
titious Heisenberg limited measurement with 10 6 atoms 
lasting one second, a "classical" projection noise limited 
apparatus would need 11 days to obtain the same level 
of precision. 

In this tutorial we attempt to provide a basic introduc- 
tion into the concept of spin squeezing. We focus on 
the theoretical concept of spin squeezing, its connection 
to entanglement and applications in quantum metrol- 
ogy. Experimental techniques to generate and detect 
spin squeezing in Bose-Einstein condensates (BEC) are 
explained in detail. Finally we shortly discuss two ex- 
emplary experiments in which spin squeezing has been 
observed in BECs. 



II. THE CONCEPT OF SPIN-SQUEEZING. 

Spin squeezing is a quantum strategy introduced in 
1993 by Kitagawa and Ueda [10] which aims to redis- 
tribute the fluctuations of two orthogonal spin directions 
among each other. In 1994 it was theoretically shown 
that spin squeezed states are useful quantum resources 
to enhance the precision of atom interferometers [11] and 
in 2001 the connection between spin squeezing and en- 
tanglement was pointed out [12] . In the field of quantum 
optics squeezed states have been proposed to overcome 
the standard quantum limit in optical interferometry al- 
ready in 1981 [HI [13]. These states are similar to spin 
squeezed states however differences arise from the dif- 
ferent variables, quadratures versus collective spin direc- 
tions [14 . The connection of spins to interferometry be- 
comes obvious when realizing that an interferometer is 
essentially a two mode system with modes \a) and \b). 
Each particle in the interferometer can be regarded as 
an elementary spin with spin a = 1/2, with two possible 
longitudinal eigenstates a z = ±1/2 corresponding to the 
two interferometer modes. The relative phase measured 
by the interferometer corresponds to the spin component 
in the plane spanned by the transversal spin directions 
<j x and <j y . The generalization of this picture to N par- 
ticles, which results in large collective spin vectors J, is 
outlined below. An extensive theoretical review of two- 
mode systems also detailing the pseudo spin description 
can be found in [T5] . 



A. Collective spins. 

The mathematical concept of a spin algebra with total 
spin J is a powerful tool to describe very different physi- 
cal systems. Any observable within such a spin J system 
can be expressed by the three spin operators J x , J y , J z 
and the identity operator. The collective spin operators 
Ji can be defined as the sum over all elementary spin 
operators (Pauli matrices) cr\ , where i = (x,y,z): 
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The three orthogonal spin components are non- 
commuting operators. Their commutation relation is 
[Ji/jj] = ieijkJk, where e^ is the Levi-Civita symbol. 
Therefore, any pair of spin operators obeys a Heisenberg 
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and A J z = ( J z ) — ( J z ) 2 is the variance in J z direction. 
A basis of the general problem can be obtained as the 
tensor product of all N bases of the individual compo- 
nents, each of dimension (2<j + 1) - the dimension of 



the Hilbert space is huge dim^A/-) = (2a + 1) N = 2 N 
and grows exponentially with the number of elementary 
spins. The length of the collective spin J is smaller or 
equal than half the number of elementary spins. The 
spin length in interferometry is usually large such that 
y/J(J + 1) « J < N/2 is a good approximation. 

In the symmetric case, that is, when all operations 
done on the spin ensemble affect each spin in the same 
way (such that there is exchange symmetry) the problem 
simplifies considerably. For example, each elementary 
spin can be prepared in the <j z — —1/2 state and maxi- 
mum collective polarization J z = —N/2 can be reached. 
Hence, the total spin length is given by J = N/2 and 
the dimension of the Hilbert space dramatically reduces 
to dim(H s ) = (2Na + 1) = (N + 1), linearly dependent 
on the number of elementary spins. One possible choice 
of a basis are the symmetric Dicke states [16] | J, m) with 
—N/2 < m < N/2. Each of the Dicke states corresponds 
to a perfectly defined particle number difference between 
the two modes \a) and \b) and - since the total number 
of particles TV is fixed - the Dicke states correspond to 
twin Fock states in the modes \a) and \b). 

Due to their exchange symmetry the elementary spins 
can be effectively described as Bosons, the Schwinger 
Bosons [17] . Employing the second quantization formal- 
ism the creation (a , b ) and annihilation operators (a, 
b) of the two modes can be related to the different spin 
components [17] 
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with n a = era and n^ = b b. Identical bosonic particles 
in two modes can be described by the symmetric spin 
model and the Schwinger representation given above can 
be used. When it comes to experiments, however, an a 
priori symmetry assumption is a strong constraint. Its 
validity has to be shown experimentally. 



1. Coherent spin states. 

Coherent spin states are the most "classical-like" pure 
quantum states of a symmetric ensemble of TV elemen- 
tary <j = 1/2 spins or of N two- mode Bosons [16] [18]. 
They are constructed by placing all N particles in the 
same single particle state, in any superposition of the 
two modes 



\o,<p) 



1 



'NV 



m{e /2)e- iv/2 a) + cos{9 /2)e iv/2 b ] ] N |vac) 



sin 



(3) 



where |vac) is the vacuum state. Any coherent spin 
state |#,(/?) can be constructed from the state |0,0) = 

-7=7 (& ) N I vac) (all N Bosons in mode \b}) by a suitable 
unitary rotation with the angles (#, tp) p~5| fTB] . 

One important aspect of the coherent spin states is 
that no quantum correlations between the particles are 
present. Therefore, a coherent spin state features equal 
variance in any direction J± orthogonal to the mean spin 
direction (#, cp) which, according to the Gaussian addi- 
tion of variances, is given by the sum of the variances 
of the 2J elementary spin 1/2 particles. The perpen- 
dicular variances A<r^ of individual elementary spins 
are by definition isotropic around (#, ip) since there are 
no subsystems that could cause any correlations (fig- 
ure [TJl) [ID]. The Heisenberg uncertainty relation ^ for 
a single elementary spin a = 1/2 pointing in a x direction 
is Aa 2 z Ady — \- \ leading to an isotropic variance of 
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for the collective coherent spin state, which identifies 

this quantum state as a minimal uncertainty state 

since (J x ) = J (see equation (I2|). We refer to the 

perpendicular spin fluctuations of a coherent spin state 
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A J_l = J/2 = N/4 as the shot noise limit. 

We go back to the first quantization formalism in or- 
der to obtain the probability distribution over different 
sets of basis states - especially the two Dicke state bases 
in the directions orthogonal to the mean spin direction. 
In order to develop a more detailed understanding of the 
coherent spin state and its fluctuations we start with the 
discussion of a special case where each particle is in a 
50/50 superposition of the two modes with 99 = rela- 
tive phase. Each elementary spin points in <j x direction 
and its quantum state is \x) = (|^, — \) + ||, +|)) /V% 
where we have chosen the Dicke states in cr z direction as 
the basis states. The probability to observe the elemen- 
tary spin in state up or down is equal |^|, ±||x)| =1/2. 
The TV atom coherent spin state is a collec- 
tion of these independent elementary spins \X) = 

[( Nh — \) + Wi + 1)) l^\ an d hence the measure- 
ment of the J z spin component is equivalent to N mea- 
surements on a single spin. The probability distribution 
over the Dicke states is binomial. We could have chosen 
equally the Dicke states in J y direction to describe the 
spin state which shows again that the spin fluctuations 
in the directions perpendicular to J x , the mean spin di- 
rection, are isotropic. 

A general coherent spin state \6,<p) in the Dicke basis 
\J,m) is described by [19 : 
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The coefficients c m (6) follow a bi- 

nomial distribution peaked around 0: 
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Due to the fundamental uncertainty in the orthogonal 
spin directions the mean direction (#, ip) of any spin state 
can not be defined with arbitrary precision. For a coher- 
ent spin state the isotropic angular uncertainty Acp = A0, 
defined by the ratio of the uncertainty of the perpendic- 
ular spin directions AJj_ to the mean spin length J, is 
given by: 
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As outlined above this limit arises as the classical sta- 
tistical limit in a system consisting of N independent 
constituents [TTJ [20] . It is this angular uncertainty which 
limits the interferometric precision in a Ramsey interfer- 
ometer. 



2. Spin squeezed states. 

Quantum correlations between the elementary spin 1/2 
particles of a collective spin J can cause anisotropic fluc- 
tuations of the collective spin in the directions perpen- 
dicular to the mean spin (figure [IJd) . In reference [10] 
quantum states are considered spin squeezed if the vari- 
ance of one spin component is smaller than the shot noise 
limit J/2 for a coherent spin state: 
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Definition (|7|) does not take the second, perpendicu- 
lar spin direction into account. Due to the Heisenberg 
uncertainty relation ^ reduction of the variance in one 
direction causes an increase of fluctuations in the other. 
Real life strategies to obtain spin squeezing might also 
involve states that are not minimal uncertainty states. 
For these states, especially for experimentally very im- 
portant non-pure quantum states, the variance in some 
other direction than the squeezed direction can be much 
larger than given by the Heisenberg uncertainty relation. 
This leads to a reduction of the effective mean spin length 

(j). 

Standard metrologic applications, for example Ramsey 
inter ferometry, require a large mean spin length. In or- 
der to measure the usefulness of spin squeezed states for 
these applications another definition of the squeezing pa- 
rameter was introduced in reference [IT] 
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whose inverse £^ measures the precision gain in a Ram- 
sey interferometric sequence relative to the standard 
quantum limit |6]). 

Spin squeezing among N constituents is related to 
many-body entanglement. In this context a third spin 
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FIG. 1. Collective spin states composed of elementary spins, a Coherent spin states. Coherent spin states in a two-level 
pseudo spin system can be thought of being composed of identical but otherwise independent elementary spins a = 1/2. 
The perpendicular variance Aa^_ — 1/4 of each elementary spin is per definition isotropic. Gaussian addition of independent 
variances result in the isotropic variance AJJ_ = N/4 for the collective spin J = N/2, b Coherent spin squeezed states. For 
coherent spin squeezed states the elementary spins have the same isotropic variance. Spin squeezing emerges due to quantum 
correlations (entanglement) between the elementary spins. The collective spin is squeezed in one direction A jj_ >min < N/4 and 
anti-squeezed in the orthogonal one A j\ max > N/4. The Q-representation is used to visualize the spins on the Bloch spheres 
(see main text). 



squeezing criterion was denned [T2] : 
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Entanglement is detected by the inequality £| < 1 as 
detailed below. Here we explicitly use the standard as- 
sumption throughout this tutorial that the mean spin 
points in J x direction and the direction of minimal vari- 
ance - if not explicitly mentioned - is the J z direction. 
£s can be used equivalently to £r to quantify spin squeez- 
ing and precision gain in interferometry and we refer to it 
as coherent number squeezing or coherent spin squeezing. 
Both quantities remain a useful parameter to detect en- 
tanglement or to measure metrologic precision gain even 
in a non-symmetric situation. To demonstrate entangle- 
ment in a collective spin system without presumptions 
not only number squeezing is required but also the co- 
herence of the quantum state has to be shown. 



3. Visualizing spin states: The Husimi Q-representation. 

Employing the Q-representation [21], collective spin 
states can be conveniently visualized on a generalized 
Bloch sphere with radius J. In order to describe the 
most general spin state, i.e. pure states and statisti- 
cal mixtures, the density matrix formalism is used [T6] . 
The density operator p in coherent spin state basis is 



given by p = f P(Q,ip)\0,(p) (Q,ip\dQ, where the inte- 
gral covers the full solid angle and dQ = sin(0) dO d<p. 
The quasi-probability distribution P{6, ip) is normalized 
to one. The Q-representation uses the diagonal elements 
of the density operator to represent the quantum state 
Q(Q, ip) = ^±1 (8,ip\p\0,ip). It gives the probability of 
finding the system specified by p in the coherent spin 
state characterized by \0,<p). The interpretation of this 
representation on generalized Bloch spheres differs from 
the usual single spin a = 1/2 Bloch sphere. In that case 
the dimension of the Hilbert space is two-dimensional 
and the quantum state representation on the surface of 
a sphere is exact. However for symmetric collective spin 
systems the dimension of the Hilbert space is 2J + 1 such 
that an exact mapping to the surface of a sphere is not 
possible. The position (0, (j>) on the spin 1/2 Bloch sphere 
describes the full quantum state, while the position on 
the generalized Bloch sphere gives only the mean spin 
direction and, within the following constraints, its fluctu- 
ations: The Q-representation projects the density matrix 
on minimal uncertainty states, in particular coherent spin 
states. The most obvious consequence is that the mini- 
mal extension of a quantum state in (0, ip) on the Bloch 
sphere is given by the uncertainties of the basis states. 
The Q-representation for visualizing collective spins is 
rather easily implemented with modern plotting tools. 
However, there exist more elaborate ways to visualize 
the quantum state of the two-mode system as described 
in [15]. The Bloch vector visualization described therein 
has the advantage that the spin noise in all three compo- 



nents of the Bloch vector is directly visible. Also the frag- 
mentation of the quantum state, that is, if the N Bosons 
do not occupy the same single particle state [15] [22] can 
be seen directly by a Bloch vector that lies inside the 
Bloch sphere. 



B. Spin squeezing and entanglement. 

Entanglement is only a meaningful concept if the sub- 
systems relative to which entanglement is defined are 
clearly specified [23-25 . Coherent spin squeezing detects 
entanglement among the elementary spins - the subsys- 
tems are the individual spin 1/2 particles [12] [26]. For 
TV distinguishable particles the definition of a separable 
state, i.e. non-entangled state, is that its iV-body den- 
sity matrix p can be written as a direct product of single 
particle density matrices pW; 
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Pk is a probability distribution to account for incoher- 
ent mixtures. Entanglement in many-body systems (for 
a general review see [23] [27] ) is defined as the non- 
separability of the density matrix p meaning the equality 
in equation (10) does not hold. 



Due to technical limitations the individual elementary 
spins can not be addressed before the measurement in 
many experiments. However, it is important to note 
that the elementary spins (the subsystems) have to be 
in principle distinguishable in order to define entangle- 
ment among them in a meaningful way [23] . For exam- 
ple, for N particles in a BEC the distinguishability is a 
priory not given. However, by a gedanken local opera- 
tion one can pinpoint each particle in space at the time of 
measurement without affecting the spin properties of the 
system. Distinguishability is now given by the position 
of each particle. If entanglement is detected in the sys- 
tem, it must have been present in the system before the 
localization, since local measurements can not generate 
entanglement [28] . Given that the atoms in the BEC are 
spaced by more than one wavelength of the detection light 
(which is usually fulfilled), this gedanken local operation 
means to overcome the technical limitations for address- 
ability and detection of the individual elementary spins. 
The fictitious pinning of the particles at the time of mea- 
surement (the photon scattering destroys the BEC) ef- 
fectively assigns a spatial mode to each particle by which 
they become distinguishable [24] [25] . The authors of ref- 
erence [29] show that the notion of entanglement in such 
a collective spin system can be generalized even without 
such a fictitious pinning of the individual particles. 
Not surprisingly different entanglement criteria can be 
found when regarding the two modes that define the 
spin 1/2, as the subsystems [30] [51] . It has been ar- 
gued that this kind of mode entanglement cannot be 
extracted in experiments without a suitable reference 
frame [32] . Recently, this problem of a missing reference 



frame for atomic fields has been addressed experimentally 
and mode entanglement has been shown in a system of 
ultracold Bosons [55] . 

Without the possibility to address the individual el- 
ementary spins before the measurement tests for entan- 
glement, that is, experimentally accessible criteria, based 
on the collective spin variables are necessary. The ob- 
servables in most experiments so far are limited to first 
and second order moments of the distributions functions 
in different spin directions due to rather small counting 
statistics and technical noise. Based on these, a com- 
plete set of inequalities, that is fulfilled for any separable 
quantum state, has been found [31] [34] [35]. Complete 
in this sense means that assuming the only information 
available are first ((J x ), (J y ), (J z )) and second moments 
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(AJ X ,AJ y ,AJ Z ) of the distribution functions. Note 
that these inequalities allow for experimental tests of en- 
tanglement, but they are no quantitative measures for 
entanglement. Such entanglement criteria only allow for 
the statement: "Yes, that state is entangled with respect 
to these subsystems." or "No, that state is not entan- 
gled with respect to these subsystems." The coherent spin 
squeezing criterion given in equation ([9| corresponds to 
a less general form of such an inequality: all separable 
states fulfill the inequality £| > 1 , but a subgroup of en- 
tangled states violate it [12] . Coherent spin squeezing as 
an entanglement criterion does not require any symmetry 
assumption. However, in the fully symmetric subspace 
many entanglement criteria simplify. For example, the 
detection of spin fluctuations in one direction below the 
shot noise limit implies entanglement [36-39 , that is 
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holds for any separable symmetric state. Equation (11) 



is identical to the spin squeezing parameter of Kitagawa 
and Ueda ^ showing that in the symmetric two-mode 
case entanglement is necessary to redistribute the fluctua- 
tions of orthogonal spin components. To avoid misunder- 
standings about the different definitions of spin squeezing 
we refer to £^ as number squeezing. 



C. Spin squeezing and interferometry. 

Entanglement in collective spin systems is not only in- 
teresting from a conceptual perspective but it has also 
been shown to provide a useful quantum resource. In 
1994 Wineland et al. [11 pointed out, that in particular 
spin squeezed states can be used to overcome the stan- 
dard quantum limit in metrology. They can be directly 
used to boost the precision of conventional Ramsey in- 
terferometers. 

Prominent applications of Ramsey interferometry [40[ 
|4T] are the definition of the time standard [5] , high preci- 
sion magnetometry [6] or inertia measurements [44] [45] . 



2T- 



")(L 



g-(Jz) 









{J.) 



input state 



after first 
n/2 pulse 



after evolution 
time r 



after second 
7r/2 pulse 
- readout - 



FIG. 2. Schematic representation of Ramsey interferometry. Panel a shows the standard Ramsey protocol depicted on the 
Bloch sphere. Beam splitters correspond to rotations of the quantum state around an axis in the equatorial plane as indicated 
by the circular shaped arrows. The sequence is described in detail in the main text. Panel b shows a similar protocol using 
spin squeezed states within the interferometer. The squeezing direction is chosen such that it finally leads to decreased noise 
in readout direction (J z direction). 



Reviews covering atom interferometry both experimen- 
tally and theoretically can be found in references [I5j l42| 



Figure [2] shows the Ramsey interferometric sequence il- 
lustrated on generalized Bloch spheres. The parameter 
S to be measured couples via a J z term in the Hamilto- 
nian resulting in a rotation of the spin around the J z axis 
during the evolution with rotation frequency S. The ac- 
cumulated relative phase p = J Sdt by which the initial 
spin state rotated around the J z direction of the Bloch 
sphere is finally converted into an observable population 
difference. 

The phase resolution of the interferometer 



Ap~ 




is determined by the projection noise AJ Z at the time 
of readout and the local slope of the Ramsey fringe 
d(J z )/dp (see figure |3]). The point of maximum sen- 
sitivity (taking finite environmental noise into account) 
is reached where the mean population difference is zero 
and the slope is maximal (d(J z )/dp) max = VN/2. The 
visibility V measures the mean spin length (J) = VN/2. 
The amount of precision gain (or loss) relative to stan- 
dard quantum limit is given by £^ , or equivalently by 

(£| ) . The measure can be expressed in visibility V 
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and spin noise in J z direction at readout A J z 
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The phase uncertainty measured as the root mean square 
deviation is Ap = ^s -jr? • 

The phase p appears in two connected contexts in this 
tutorial. On the one hand, for coherent spin states, it 
is a variable characterizing the full quantum state of the 
system (together with 0). On the other hand, it has the 
meaning of a relative phase between the two modes of 
the interferometer (or more general between the modes 
| a) and \b) of the two-mode system). In both cases the 
phase p can be interpreted as the mean angle of the 
spin state in longitudinal direction (that is, in the J X1 
J y plane) on the Bloch sphere as schematically shown in 
figure |2j The angular uncertainty in this direction finally 
limiting th e precis ion of a Ramsey interferometers is Ap 
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Spin squeezed states feature reduced noise in one of 
the spin directions but excess noise in another direction 
can be present either due to a non-Heisenberg limited 
quantum state or due to an incoherent mixture of several 
quantum states. The former might limit precision in 
standard Ramsey interferometry, but specific correlated 
quantum states enable even enhanced interferometric 
precision in a generalized interferometer [48]. The latter 
is easily limiting interferometric precision at a level 
above the standard quantum limit and experimentally 
it requires a large effort to prevent decoherence due to 
technical noise from the environment or due to finite 



1 For two-mode systems a third notion of a relative phase via the 
definition of phase eigenstates exists that will not be used in this 
tutorial [151 [46)1471 . 




FIG. 3. Precision limit in Ramsey interferometry. The phase estimation precision in Ramsey interferometry using a coherent 
spin state (gray) and a spin squeezed state (black) are schematically explained. The main figure shows a Ramsey fringe whose 
visibility V is maximal for a coherent spin state (V = 1) but smaller for a spin squeezed state. Nevertheless the phase precision 
for a squeezed state outperforms the precision obtained for classical interferometers as shown in the zoom around the point 
of highest sensitivity. The projection noise is suppressed for the spin squeezed state such that the ratio of projection noise 
and slope of the Ramsey fringe is smaller by a factor ^s compared to the standard quantum limit, which explains the gain in 
inter ferometric precision. 



temperature in the system. Large noise - quantum or 
classical - even in a spin direction that is not directly 
measured has a degrading effect on interferometric 
precision in a standard interferometer which arises due 
to the curved surface of the Bloch sphere. As soon as 
the noise amplitude is large enough such that the area of 
uncertainty can no longer be approximated by a plane, 
the mean spin is effectively shortened and the visibility 
decreases V < 1. 



Entanglement can be used as a quantum resource in 
a Ramsey interferometric sequence in different ways. In 
order to increase the phase sensitivity either the slope of 
the signal d(J z )/dcp has to be increased or the projec- 
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tion noise A J z has to be decreased. Slope increase can 
be reached by Schrodinger cat type entanglement which 
involves maximally entangled states that are very fragile 
to decoherence. Therefore they have been realized so far 
with very few particles only [49H5T] . 
Spin squeezing aims to decrease the projection noise. 
This is possible in gradual steps meaning that depending 
on the amount of spin squeezing the precision is gradu- 
ally increased. Therefore - at least for moderate levels 
of spin squeezing - these states are less fragile and they 
have been realized with a large number of particles but 
only with a relatively small squeezing factor [4, 52-60 . 
Interferometric sensitivity for a coherent spin state and 
a coherent spin squeezed state is illustrated in figure [3] 
For the coherent spin squeezed state the decreased quan- 
tum fluctuations AJ Z reduce the projection noise while 
the increased fluctuations A J y cause a slight decrease of 



the mean spin length and therefore of the visibility of 
the Ramsey fringe. Nevertheless, the ratio of projection 
noise and slope of the Ramsey fringe - and therefore the 
phase sensitivity - is increased. 



D. Heisenberg limit in quantum metrology. 

The ultimate limit for metrologic precision is the 
Heisenberg limit [7 j, where the phase estimation error 
Aip is given by Aip = -^ for N resources used in a sin- 
gle measurement. This fundamental limit can in princi- 
ple be reached with both approaches mentioned above, 
Schrodinger cat type entanglement or spin squeezing. 

In the context of quantum metrology the Schrodinger 
cat state is sometimes called a NOON state [61 . Its name 
originates from its form in Fock states basis \NOON) 



(\N,0) 



\0,N))/y/2. It is a coherent superpo- 



sition of all atoms in state \a) and zero atoms in state 
\b) and vice versa. In spin representation the NOON 
state is the superposition of the two maximal Dicke states 
\NOON) = (| J, -J) + e^" | J, J))/y/2. The increase of 
the signal slope for a NOON state is obvious since the 
phase acquired between the two components ip^ = Nip is 
N times larger than for a coherent spin state [20, 62] [63]. 
Experimentally it is important to note that the readout 
of the interferometer can not be realized by measuring 
(J z )- The reason is the vanishing mean spin length (J x ) 
of this state. It has been shown that the parity of the 
state is a useful experimental observable to make use of 
NOON states in interferometry and to reach the Heisen- 



berg limit [491164]. 

Further proposals to achieve Heisenberg limited inter- 
ferometry involve twin-Fock states \N/2,N/2) [65], [66]. 
The mean spin length for these states is also vanish- 
ing such that special readout schemes are required [67] . 
Recently such states have been produced experimen- 
tally [33] [68] [69] and their usefulness for quantum en- 
hanced metrology has been shown [69] . 

Spin squeezed states allow for the optimization of in- 
terferometry gain demanding a finite mean spin length 
such that standard interferometric readout can be used. 
The optimum £r for a given mean spin length was found 
numerically in reference [26]. An experimental protocol 
to generate spin squeezed states close to the Heisenberg 
limit was proposed in reference [70j [71]. The authors 
show that the Heisenberg limit can be reached asymp- 
totically in the number of atoms. 

The precision of Ramsey type interferometric measure- 
ments increases with evolution time within the interfer- 
ometer (cf. section [I]). The same inter atomic interac- 
tions, that can be used to generate entangled states ini- 
tially and to overcome the standard quantum limit (as 
shown in this tutorial), however, might spoil the preci- 
sion of the interferometer if present during this evolution 
time [72] . Hence, interaction control (for example us- 
ing Feshbach resonances) is crucial for future applications 
that use standard interferometric readout. 

Recently it has been emphasized that the Fisher 
information is the most general criterion to measure 
phase sensitivity [48, 73 since it saturates the quantum 
Cramer- Rao bound [2j [74] [75]. Calculating the Fisher 
information for a coherent spin state state evolving un- 

~ 2 

der the non-linear Hamiltonian H = \Jz > where \ 
parametrizes the nonlinearity, Heisenberg like scaling for 
the phase precision was found [48] after long evolutions 
times where the potential precision gain in a standard in- 
terferometer was already decreased. The quantum state 
here is neither necessarily a NOON state nor a coherently 
spin squeezed state. However standard interferometric 
readout can not be used to extract the phase information 
and a new type of Bayesian readout might be employed 
which was experimentally demonstrated in [75] . 



III. PREPARATION OF ATOMIC SPIN 
SQUEEZED STATES. 

Strong effort has been put into the experimental gen- 
eration and detection of spin squeezed states in large col- 
lective spin systems. Early experiments detected spin 
noise beyond the shot noise limit in a cold ensemble of 
Cesium atoms prepared by quantum state transfer from 
non-classical light to the atoms [76], or in a Cesium va- 
por cell at room temperature by quantum non-demolition 
measurement [77]. Indirect measurements of number 
squeezing of ultracold atoms in optical lattices and dou- 
ble well systems followed, however, the direct detection of 
the fluctuations was not possible [78-83 . Several exper- 



iments generating and detecting coherent spin squeezing 
of large collective spins in atomic gases were done in last 
four years [53] [Ml 133 W\ [59] [60] [84] . Full Ramsey inter- 
ferometers working beyond the standard quantum limit 
have been realized using coherent spin squeezed states 
within the interferometer [55, 60] 155]. 



A. Spin squeezing with Bose-Einstein condensates. 

Different methods have proven useful to generate spin 
squeezed states in atomic samples. Many of them in- 
volve light, for example by quantum state transfer, quan- 
tum non-demolition measurements or cavity induced in- 
teractions [86] . In contrast, ultracold quantum gases of- 
fer the possibility to directly use the intrinsic interaction 
between the atoms. These interactions lead to nonlin- 
ear terms in the Hamiltonian necessary to generate spin 
squeezed states [12j [87] . 

We will take intrinsic s-wave interactions among atoms 
in a two-mode BEC as an example here and discuss how 
spin squeezing can be generated. Examples for such two- 
mode systems based on internal or external degrees of 
freedom are shown in figure [4j In such two- mode BECs 
the inter atomic interactions lead to a nonlinear term in 

A 2 

the Hamiltonian proportional to J z [22, 88, 89 . A re- 
view on the experimental realization of the external two 
mode system can be found in [46] and for the internal 
system details are given in [90] [91] . The theoretical de- 
scription for both cases, including the possibility for frag- 
mentation of the quantum state [92] , has been reviewed 
in detail recently [15 J. Thus, we only sketch the con- 
nection between the general Hamiltonian and the pseudo 
spin description for the external double well case in the 
following. The general Hamiltonian is 



H 



I 



2m 



drdr / xj/ t ^ /t V(r-r / )M / 



(14) 



from which, by making the two- mode ansatz \I/(r) = 
a (r)a + (\)\) (r) 6, the two-mode Bose- Hubbard Hamilto- 
nian is obtained (for improved readability we did not 
write the spatial dependence of the operators in equa- 
tion (14) explicitly). 



H, 



BH 



hx 



(tfa-tfb) 2 



Here we defined three parameters in the Hamiltonian, 
the nonlinearity x characterizing the inter atomic inter- 
actions, the Rabi type coupling between the two modes 
Q and a differential energy shift S. The remaining step is 
to write the Hamiltonian (15) in spin language by using 
the Schwinger representation introduced in equation (|3| : 



H = —h£iJ x + h\J z — hSJ z 



(16) 




FIG. 4. Examples for two- mode BECs. a A BEC in a double 
well potential for sufficiently high barrier between the two 
wells can be approximated as a two-mode system where the 
mean field modes in the two wells (j) a and cj>b define the modes. 
b Another example for a two-mode BEC is a spinor BEC with 
two populated hyperfine spin states \a) and \b) and identical 
spatial wave function 0. 



Note that the Hamiltonian in both cases, external and 
internal degrees of freedom, can be written in the form of 
equation (fl6|) . The actual calculation of the parameters 



involves spatial integrals over the wave functions and 
differs in both cases. In this context it is important to 
note that the spatial wave functions <p are not the single 
particle ground state wave functions in the confining 
potential, but the mean field wave functions obtained by 
coupled Gross-Pitaevskii equations for the two modes 
(For several different approaches see [TSJ [46j [92] [93]). 
Details on the calculation of the parameters in the 
Hamiltonian can be found for example in [90, 91, 94] [95]. 

Three regimes can be distinguished by looking at the 
quantum fluctuations in the orthogonal spin directions 
for the ground state of the Hamiltonian in equation ( 16 ) 
(see figure [5|. For the following discussion we assume 
an unbiased situation, that is 5 = 0. The important 
parameter measuring the relative importance of tunnel 
coupling and interaction is A = Nx/Q- This can be 
seen directly by normalizing the spin operators to spin 
1/2 such that the total number of atoms appears explic- 
itly. Then A is the ratio of the nonlinear to the lin- 
ear coefficient. The Rabi regime is the noninter acting 
limit where A <C 1. Here the ground state is a coherent 
spin state with maximal mean spin length ( J x ) — N/2 
and equal spin fluctuations in the orthogonal directions 

~ 2 - 2 

AJ y = AJ Z = N/4. For increasing interactions fluc- 
tuations in the relative atom number (equivalent to fluc- 
tuations in J z ) become energetically unfavorable and 
the system enters the Josephson regime (1 < A < N). 
Its name originates from the analogy of the Hamilto- 
nian ( [16] ) to the Josephson Hamiltonian in solid state 
physics|22j [96] [97] . The ground state is a coherent spin 
squeezed state with reduced spin fluctuations in J z direc- 
tion at the cost of increased fluctuations in J y and hence 
a reduced mean spin length. The tunnel coupling is neg- 
ligible in the Fock regime (A >> N). The ground state 
is an extremely spin squeezed state with vanishing mean 
spin length. It finally approaches the twin Fock state, an 



eigenstate of J z with eigenvalue zero. 

There are two possibilities to control the regime pa- 
rameter A. For a given total atom number N either 
the tunnel coupling Q or the nonlinearity \ nas to be 
changed. While the former is straight forward, there 
are again several possibilities to change the nonlinearity 
X- The rather direct approach is to control the scatter- 
ing properties [60 , which is directly possible when using 
quantum states with accessible Feshbach resonances [98] . 
The second method, which has also been realized experi- 
mentally, is to control the wave function overlap between 
the two states [841193] . 

Assuming the weakly interacting regime can be real- 
ized experimentally there are two strategies to generate 
coherent spin squeezed states. \^j The first is an adiabatic 
approach where the regime parameter is changed slowly 
driving the system deeper into the Josephson regime [70] . 
For some experimental settings a fast non-adiabatic ap- 
proach might be preferable [10] [12] where spin squeezing 
develops during the evolution after a quench in A. In the 
following we discuss the advantages and disadvantages of 
the two strategies. 



1. The adiabatic approach. 

One of the advantages of ultracold atom systems is the 
ability to control the parameters of the Hamiltonian dy- 
namically. Starting from a strong tunnel coupling the 
regime parameter A can be adiabatically changed. Such 
an adiabatic transformation of the quantum state is a ro- 
bust way to generate coherent spin squeezed states. The 
tunnel coupling is present throughout the transforma- 
tion, making the system less prone to phase noise (fluc- 
tuations in 5). This can be intuitively understood by the 
concept of the spin echo technique [100 : For small differ- 
ential energy shifts 5 <C ft the continuous tunnel coupling 
causes a continuous sequence of ir echo pulses counteract- 
ing the dephasing. Another advantage of the adiabatic 
technique is its usefulness in a finite temperature environ- 
ment, where the initial state is a thermal mixture of the 



lowest many-body eigenstates of the Hamiltonian (16). 



In the Josephson regime the adiabatic increase of A re- 
sults in an effective "cooling" of the system [59] . The term 
cooling should be used with care here since no entropy 
is removed from the system. With increasing parameter 
A each of the lower eigenstates of the system features re- 
duced fluctuations in J z , however at the cost of increased 
fluctuations in J y and a loss in coherence. Nevertheless, 
the gain in J z direction outweighs the loss of coherence 
for small enough temperatures, such that the coherent 
spin squeezed regime can be reached by this method. 



2 Later in this manuscript we will discuss the effects of finite tem- 
perature on spin squeezing in a double well setup, where the 
advantages of starting in a weakly interacting regime are dis- 
cussed. 
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FIG. 5. Ground state spin fluctuations. Panel a, shows the representation of the ground state of the collective spin Hamil- 
tonian (16) on the Bloch sphere in the three different regimes. The histograms represent the probability distribution of the 



quantum state over the J z and J y Dicke eigenbasis respectively. The distribution in J z narrows while the variance in J y in- 
creases. Note that the distribution of J y in the Fock regime shows fringes with a period of 1/N which in principle can be used 
for increased sensitivity in interferometry [481 199) . b, shows the coherent spin squeezing parameter £| as a function of A. The 
calculation was done for 100 atoms. 



It has been theoretically shown that the Heisenberg 
limit can be reached asymptotically in the number of 
particles by the adiabatic technique [70j [71] . Experimen- 
tally, however, it is a major problem to assure adiabatic- 
ity during the full range of the parameter change. At 
least when changing the regime parameter by tuning the 
tunneling coupling (for fixed and typically small nonlin- 
earity x) the required timescale for adiabaticity eventu- 
ally diverges setting a limit to the attainable coherent 
spin squeezing factor £# • The long times required to gen- 
erate coherent spin squeezing adiabatically are the main 
disadvantages of this approach. 



2. The diabatic approach. 

For some experimental systems the robust but slow 
adiabatic approach to generate spin squeezing might 
not be optimal. A diabatic technique, in which the 
Hamiltonian parameters are quenched (changed non- 
adiabatically) , can produce spin squeezing in a much 
shorter time. The idea of the diabatic one-axis twisting 
was introduced in 1993 by Kitagawa and Ueda. Starting 
with a coherent spin state with mean spin in J x direction 
the tunnel coupling is suddenly switched off such that the 
nonlinear term in the Hamiltonian governs the evolution 



U{t) = e~ itxfz [HHH21I93] . The nonlinearity causes a 
shearing of the uncertainty region on the Bloch sphere 
since the rotation speed around J z is itself J z dependent. 
The scheme is similar to Kerr effect based squeezing pro- 
tocols in quantum optics where a material with a finite 
Kerr nonlinearity is used such that the refractive index 
flight oc n2 |£^| is proportional to the light intensity \E\ . 
The light experiences intensity dependent phase modula- 
tion within this medium resulting in quadrature squeez- 
ing [TOT] . 



The one axis twisting evolution of an initial coher- 
ent spin state can be nicely visualized. Earlier in this 
manuscript it has been discussed how a coherent spin 
state on the equator of the Bloch sphere \0 = 7r/2, ip = 0) 
can be described as a coherent superposition of several 
Dicke states where the probability distribution over these 
basis states is binomial. Each Dicke state |J, m) com- 
posing the coherent spin state rotates with a different 
frequency around the J z axis where the difference in ro- 
tation frequency between next neighboring Dicke states 
is x (figure [6k). For short evolution times r < Td ep h 
this shearing effect results in coherent spin squeezing un- 
der an axis rotated by the angle ao(r) with respect to the 
equator of the Bloch sphere (figure pp). Td ep h = (^mX) -1 
is the dephasing time, after which the mean spin length 
2(J)/N has dropped from a value close to unity to e _1 
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FIG. 6. Non-adiabatic spin squeezing by one axis twisting. Panel a illustrates the initial coherent spin state on the equator 
of the Bloch sphere. The histogram shows the binomial probability distribution over the Dicke states. The zoom window 
illustrates the quantization on the J z axis - the different Dicke states - of which each eigenstate m rotates with a different 
angular frequency mx around the vertical axis. Panel b illustrates the quantum state after a short evolution time r < Tdeph- 
The isotropic uncertainty developed into an elliptical one with spin squeezing present under an angle ao . The histogram shows 
the squeezed probability distribution over the eigenstates in a coordinate system rotated by ao. 



due to the interaction induced spread of the state around 
the Bloch sphere [102 . In general Td ep h is inversely pro- 
portional to the extension of the quantum state over the 
Dicke basis <j m which, for a coherent spin state on the 
equator of the Bloch sphere, is a m = y/N/2. 
At evolution times r > Td ep h reduced quantum fluctua- 
tions under a certain axis are still present but the mean 
spin length is strongly reduced such that the quantum 
state is no longer optimally coherently spin squeezed. 
After even longer times the dynamics show highly non- 
classical interference effects [103-105 finally resulting 
in a revival of the mean spin length at r rev = 2nx~ 1 
when each neighboring pair of Dicke states is in phase 
again |104j . The best achievable noise suppression in 
general interferometry still increases with evolution time 
even if coherent spin squeezing degrades. It has been 
recently shown [48] that a new type of Bayesian inter- 
ferometer readout can be employed to make use of these 
quantum states. As typical for macroscopic superposi- 
tion states the highly entangled states occurring for long 
evolution times are extremely sensitive to decoherence. 
For example, particle loss leads to an exponential decay 
of the revival contrast with a N times enhanced decay 
constant [106 . 

The maximum squeezing achievable with the diabatic 
one axis twisting technique is £| « TV -2 / 3 [T0j [12]- This 
has to be contrasted to the maximum squeezing that can 
be generated using the adiabatic technique which ap- 
proaches the Heisenberg limit £| ~ iV _1 at the bound- 
ary to the Fock regime where the mean spin length is still 
reasonable high [70]. The maximum squeezing is better 
in the adiabatic case however the adiabaticity criterion 
requires long evolution times r a di a b ~ 1/x- The best 



squeezing in the one axis twisting protocol is achieved 
after r^ a « 7V -2//3 /x, a factor of 7V 2//3 faster favoring 
the diabatic protocol for certain applications. 



The generation of spin squeezing might be optimized 
using different Hamiltonians. Stronger maximal squeez- 
ing, for example, can be reached by a two axis twisting 
scheme [10]. Another example, optimizing the required 
time to achieve squeezing, is to use evolution under the 
Hamiltonian ( 16 ) after a quench to a small but finite tun- 



nel coupling (Into the Josephson regime) [107] . For ini- 
tial states with relative phase (p = ir the evolution leads 
first to coherent spin squeezing and for longer times to 
Schrodinger cat like states [108] 1109] . 



When discussing the maximum achievable coherent 
spin squeezing one has to keep in mind that with de- 
creasing £| the generated states become more and more 
sensitive to decoherence. For the exemplary case of the 
one axis twisting evolution this problem has been stud- 
ied in detail taking into account particle loss, random de- 
phasing (for example due to noise in £, see equation ( 16)) 
and multi mode effects due to finite temperature which 
are always present in real condensates (11014112] . Includ- 
ing these effects the achievable coherent spin squeezing 
is no longer proportional to TV -2 / 3 but tends to a fi- 
nite but small constant value in the thermodynamic limit 
N —> oo, showing that a large metrology gain can still be 
reached [111] . Reference jl 12] gives a review on decoher- 
ence effects in dynamically generated spin squeezing with 
BECs and points out the effective dephasing effect of the 
different sources of decoherence. 
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FIG. 7. Ramsey fringe using a coherent spin squeezed state. The figure shows a Ramsey fringe obtained for a coherent spin 
squeezed state prepared using the one axis twisting sequence for an evolution time r <C Td ep h- The solid line is a sinusoidal fit 
through the data from which the visibility (here V = 1.00 ± 0.02) and mean spin length (J) = NV/2 can be deduced. This 
figure was originally published in the supplementary online material to [60] . 



IV. DIRECT OBSERVATION OF COHERENT 

SPIN SQUEEZING - EXPERIMENTAL 

CHALLENGES. 

Both approaches the adiabatic and diabatic one out- 
lined in the last section require precise control over the 
experimental parameters. Especially so far neglected un- 
controlled differential energy shifts 5 are the main exper- 
imental challenge to observe spin squeezing. Another im- 
portant issue, mainly in experiments with BECs in dou- 
ble well traps, is finite temperature which easily limits 
the observable spin squeezing. The detection of atomic 
noise at the shot noise level in a small optically thick 
cloud is challenging and requires careful calibration of 
the imaging system. Furthermore, a method to detect 
the mean spin length has to be available experimentally 
in order to detect coherent spin squeezing. In this sec- 
tion we will discuss these challenges in detail, starting 
with the detection issues. 



A. Measurement of the mean spin length. 

When working with internal degrees of freedom the 
measurement of the mean spin length can be done by a 
Ramsey experiment [40j [41] . A Ramsey scheme is usu- 
ally relatively easy to implement given the availability of 
precise microwave and radio frequency sources (or Raman 



lasers) to control the tunnel coupling |100j . This coupling 
Q usually can be made much stronger than the nonlin- 
earity \ sucn that the Rabi regime is reached. When 
working with two photon transitions, which is often re- 
quired to use states less sensitive to magnetic field fluctu- 
ations, care has to be taken to accurately determine the 
resonance. The resonance is shifted during the coupling 
with respect to the free evolution times due to the ac- 
Stark shift induced by the coupling fields [90] [91]. The 
readout of the Ramsey interferometer requires precise de- 
tection of the atom number difference in the two internal 
states over the full dynamic range. The mean spin length 
(J) = NV/2 is then inferred from the visibility V of the 
Ramsey signal. An example for such a Ramsey fringe 
revealing a visibility of V = 1.00 ± 0.02 is shown in fig- 
ure [7] [60]. 

For external degrees of freedom, for example in a dou- 
ble well trap for BECs the mean spin length can be re- 
vealed from the interference pattern after time of flight 
expansion [59, 96] I113( II 14] . Additional complications 
might arise from a finite imaging resolution or due to 
excitations of additional spatial modes when expanding 
the condensate, both of which limit the visibility of the 
interference pattern in a single shot. Care has to be 
taken when extracting the average visibility V. In gen- 
eral the analysis method chosen has to include the effects 
of the possibly reduced single shot visibility, which might 
also be caused by multi mode effects in the initial dou- 
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FIG. 8. Calibration of the imaging system, a The correction factor c gpe is determined by best matching the observed profile to 
the numerical simulations. The figure shows an exemplary calibration, b The Clebsch-Gordan coefficient is obtained in a second 
step by requiring imaging independent atom number measurement. The figure shows an example of a successful calibration. 



ble well [90, 91 . Exemplary interference measurements 
for a BEC released from a double well trap are shown in 
figure [Tffi later in this manuscript. 



B. Absolute atom number measurements at the 
shot noise level. 

The detection of number squeezing requires an accu- 
rate calibration of the deduced total atom number and 
the linearity of the imaging system. Simultaneously a 
high spatial resolution is often required [59j [60] . To re- 
alize sufficiently large nonlinearities x m the mesoscopic 
atom number regime a high atomic density and there- 
fore tight confinement of the atoms is necessary, lead- 
ing to optically think samples. This complicates the use 
of in principle single atom sensitive fluorescence detec- 
tion methods |115[ 1116] which suffer from light induced 
collision losses if the individual atoms are not spatially 
separated during the detection. An alternative is to use 
absorption imaging [1 13] in the high saturation intensity 
regime [11 7] which provides good signal to noise ratio 
and high spatial resolution, however, its calibration is 
non trivial as outlined below. 



1. Calibration of the imaging system. 

In high intensity absorption imaging [11 7] the full Beer- 
Lambert absorption formula is necessary to calculate the 
atomic column density riij from the counts per pixel 



7>P!C,i,J 



on the picture containing the absorption infor- 
mation and on the reference picture iV^^y (i,j indexes 
the CCD camera pixels but it is mostly omitted below 
to ensure better readability). For the image analysis the 
mean light intensity on the reference picture has to be 



normalized to the mean light intensity on the picture 
with the atomic signal. For its determination the area 
containing the atomic signal is not taken into account. 
Not only the optical density Od — li^-W^ref/A^pic) but 



also the difference in the counts A = (N, 
contribute to the signal: 



7,ref 



N, 



7, pic 



d 1 l n 

n= — ( — O d 

&0 Ccg 



Cccd C| 



gpe 



A) 



(17) 



Here d is the linear extension of the CCD pixel taking 
the magnification into account, do = 3X 2 /2tt the resonant 
cross section, c cg the Clebsch-Gordan coefficient, c gpe a 
correction factor obtained from a comparison to simula- 
tions as explained below and r is the imaging pulse dura- 
tion. The factor c cc d = hu)/(d 2 rjQI sa t) contributes to the 
linear part of the formula where Q is the total quantum 
efficiency and 77 the gain factor of the camera, uj is the 
light angular frequency and I sat the saturation energy of 
the transition. The gain factor 77 is determined by mea- 
suring the noise curve AiV 7 ((iV 7 )) of the CCD camera. 
It can be measured over the full dynamic range using an 
incoherent light source (e.g. a LED) such that the light 
hitting the active region of the CCD is uncorrelated fea- 
turing shot noise limited noise characteristics. The factor 
r\ follows from the slope of the CCD's noise curve via a 
linear fit in the working region where the camera noise is 
not read out noise dominated. 

The calibration of the c cc d factor is based on the quan- 
tum efficiency measurement which itself relies on the 
power meter calibration. Furthermore, for a typical nu- 
merical aperture around 0.5 approximately 5% of the 
solid angle is covered by the imaging lens such that a 
non negligible fraction of the scattered photons reaches 
the camera. This effect is not taken into account in equa- 
tion (17). A possible method to cross check the inferred 



number of atoms is to prepare the BEC in a very well 
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FIG. 9. Nonlinear effects in absorption imaging of small 
atomic clouds. For small clouds (gray boxes) and up to in- 
termediate imaging intensities (in the optimal signal to noise 
regime) the atom number is systematically underestimated 
with rising optical density. Linearizing the result around 750 
atoms yields a factor of two smaller slope while the absolute 
number is underestimated by only 15%. This results in a 
strong bias of the measured atom number fluctuations in this 
regime. The reference atom number in this figure was mea- 
sured for 400 fis expansion time and with high imaging inten- 
sity ///sat (black circles) [59, 90, 91 . This figure was origi- 
nally published in the supplementary online material to 59 . 



known trap configuration and image it with very high 
imaging intensity (typically 50/ sat ). In this regime only 
the linear part of the Beer-Lambert formula contributes, 
but the signal to noise ratio is poor. Imaging is done with 
a very short pulse (2 /is) to minimize diffusive broaden- 
ing of the profile due to photon scattering. The extracted 
average profile along an axis much bigger than the opti- 
cal resolution is calculated for different correction fac- 
tors Cgp e and compared to theoretical profiles obtained 
from three dimensional Gross-Pitaevskii equation simu- 
lations varying the total numbers of atoms (figure [8k). 
The correction factor is finally found by minimizing the 
quadratic deviation between the simulated and measured 
profile. The expected corrections are small such that c gpe 
is should be close to unity. 

Calibration of the imaging is completed by the mea- 
surement of the Clebsch-Gordan coefficient c cg . This is 
done by imaging the atoms with varying imaging inten- 
sities I/I sat . A correct Clebsch-Gordan factor means in- 
tensity independent atom number measurements (see fig- 
ure |]b). 

The total atom number follows from direct summation 
over the region where the atoms are detected (not by 
a fit which might bias fluctuation measurements). The 
size of this region is typically chosen to three standard 
deviations as obtained from a Gaussian fit to the atomic 
cloud. 
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FIG. 10. Correctly calibrated imaging system. The vari- 
ance of a coherent spin state on the equator of the Bloch 
sphere is measured versus the total atom number to test the 
imaging calibration. The linearity (with a fitted slope close 
unity, black line) of the measured data confirms the indepen- 
dently obtained calibration of the imaging system. This figure 
was originally published in the supplementary online material 
to EDI. 



2. Non-linear effects. 

Quantitative measurements of the atom number using 
absorption imaging is limited to atomic clouds with a size 
larger than the optical resolution of the imaging system. 
For to small clouds the non-linear part of equation (fl7|) 



causes a systematic error in the detected atom number. 
This can be easily understood in the limit where the op- 
tical resolution is much larger than the pixel size. The 
detected column density n per pixel is calculated from 
photon counts integrated over the size of the pixel 



oc In 



7. TV. 



7,ref 



Jpix ^7,Pi< 



(18) 



where we omitted the linear term in equation (17). Cor- 
rect atom number calculation n true however requires to 
calculate 



n tr 



oc 



^pix \^7,pic/ 



(19) 



In general n ^ n true holds such that the calculated atom 
number is wrong. When the pixel size is smaller than the 
point spread function / of the imaging system the same 
argument holds, but the main averaging effect is due to 
the convolution with the point spread function, e.g. the 
replacement iV 7 — )> (7V 7 */) has to be made, where (g* f) 
means the convolution of the functions / and g. Figure [9] 
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FIG. 11. Spectral sensitivity to phase noise. The figure shows the normalized spectral phase noise sensitivity of a Ramsey 
sequence for an exemplary free evolution time of 20 ms (gray line) . The black line shows the altered sensitivity including a spin 
echo after 10ms. 



shows an example of this effect. The nonlinearity in the 
atom number calculation for high optical densities and 
small atomic clouds is clearly visible. It leads not only 
to an underestimation of the absolute atom number but 
also to an underestimation of the fluctuations in the high 
density region due to the smaller local slope. By slightly 
expanding the atomic clouds before imaging this effect 
can be minimized. 



3. Photon shot noise. 

The atoms are detected via their resonant interaction 
with the probe light whose noise characteristics add to 
the atomic noise of interest. At least the shot noise of the 
detection light adds to the detected noise, but especially 
in absorption imaging interference fringes easily increase 
this minimal background noise level. These fringes orig- 
inate from motion of the imaging systems optics or air 
motion in the light path. Special care has to be taken 
to eliminate these effects. The detection light shot noise 
contribution can be inferred pixel by pixel on both, the 
reference and absorption picture, using the camera noise 
calibration curve. Standard error propagation of equa- 
tion (17) allows for a conversion of the photonic variance 
AiV 2 ™ V /w^ i n to atomic variance: 



AiV 2 • • 

psn,i,j 




Here the photon shot noise contribution AiV 2 i • is cal- 
culated by a expansion up to second order, assuming 
a Gaussian distribution for the photon statistics p{N 1 ). 
Photon shot noise on different pixels (i,j) is uncorrelated 
such that the total contribution in a given area of the 
CCD chip AiV 2 sn can be obtained by summation over the 
variance per pixel. The photon shot noise is independent 
of the atomic noise AiV 2 . Thus, the best estimate of the 
true atomic fluctuations can be obtained by subtracting it 
from the detected fluctuations AiV 2 = AiV| et - AiV 2 sn . 
The uncertainty in the estimated amount of photon shot 
noise is typically in the order of a few percent. 
In the strong saturation regime the imaging light inten- 
sity / controls the transparency of the atomic cloud [117j . 
The imaging signal to noise ratio can be optimized by 
choosing the proper light intensity depending on the 
atomic density and the camera parameters. 



4- Independent tests of the imaging calibration. 

From the discussion above it is clear that correct cal- 
ibration of the imaging system is non-trivial. Therefore 
it is important to have an independent test of the de- 
duced atomic variances. One possibility is to hold an 
atomic sample in the trap until a large fraction of the 
atoms is lost by background gas collisions. This uncor- 
related one-body loss tends to restore a Poissonian atom 
number distribution. However, the Poissonian limit is 
only reached asymptotically and the fluctuations might 
be super-Poissonian for an initial distribution governed 
by experimental instabilities or slightly sub-Poissonian 
if the loss is dominated by two or three body colli- 
sions [59J [901 El EH EH] • 
When working with internal degrees of freedom coherent 
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FIG. 12. Position noise leads to relative atom number fluc- 
tuations. The schematic illustrates the connection of number 
difference fluctuations and relative position fluctuations Az 
between optical dipole trap and optical lattice. In equilibrium 
the chemical potential \i is the same in both wells which re- 
sults in an occupation number difference between the left and 
right mode dependent on the differential energy shift AVharm 
between them. 



spin states can be prepared by precise electromagnetic 
coupling pulses. Their well known noise characteristics, 
that is, relative atom number fluctuations at the shot 
noise limit, can be directly used for the imaging calibra- 
tion. Figure [To] shows such a calibration check which 
reveals a correctly calibrated imaging over a large range 
of total atom numbers [60, 9Q]l9T]. 



C. Differential energy shifts. 

In the discussion of the squeezing generation we as- 
sumed vanishing external noise during the evolution. In 
real experiments, however, fluctuations of the differential 
energy shift S coupling via the J z term in the Hamilto- 
nian (16) cause phase noise [112 . Both, the adiabatic 



and non-adiabatic approach to generate coherent spin 
squeezed states suffer from this noise problem. Especially 
the non-adiabatic method is very sensitive to DC phase 
noise (fluctuations from shot to shot). A spin echo tech- 
nique |100j can be used to effectively remove the sensitiv- 
ity to these shot to shot fluctuations, however, the sen- 
sitivity to low but finite frequency fluctuations remains 
problematic. In general the spectral sensitivity is depen- 
dent on the duration of the free evolution and knowing 
the noise spectrum a specially designed spin echo tech- 
nique can be used to optimize the spectral sensitivity. 
Figure [TT] shows an exemplary sensitivity curve versus 
frequency of the environmental noise with and without 
spin echo. 

In the adiabatic case the final relative atom number 
fluctuations are also influenced by the noise in the differ- 
ential energy shift, however, the main contribution is due 
to shot to shot fluctuations. The continuous coupling of 
the two states results effectively in continuous spin echo 
pulses counteracting dephasing. The remaining DC sen- 
sitivity can be understood intuitively for the example of 



the external double well system, but the argumentation 
also holds for an analogous internal system. The sys- 
tem response to an imbalance in the offset energy for the 
two wells (a DC differential shift) is to adapt the relative 
atomic population such that the overall chemical poten- 
tial is constant. This means an increased population of 
the lower state such that the extra interaction energy bal- 
ances the external energy offset. An estimation for the 
noise sensitivity can be obtained by the following con- 
sideration where the main experimental challenge, the 
relative movement of the dipole traps that generate the 
trapping potential j46[ I12QJ becomes obvious. Figure 12 
illustrates this situation. 

The optical dipole trap generates harmonic confine- 
ment in axial direction Vharm = muj^z 2 /2, where m 
is the atomic mass. Fluctuations of the energy differ- 
ence AVharm between the two potential minima sepa- 
rated by d due to fluctuations of the barrier position 
Az is given by AVharm = dv £ a z Tm Az = muj 2 z dAz. The 
overall chemical potential in local density approxima- 
tion [121 is \i = Vharm(^) + /J>at(z) an d in equilibrium 
it is the same for the two wells. The contribution due 
to inter atomic interactions fJL a t( z ) balances the change 
in VharmO): A/i a £ = ^An a = 2xAJ 2 , where we used 

-q^- — 2x [SS] an d J z = (rib — n a)/^- To give an es- 
timation of the required position stability of the trap- 
ping potentials we calculate the actual numbers for the 
trapping parameters of the double well trap in refer- 
ence [59 . The nonlinearity is 2\ ~ 27r • 1 Hz and the 
longitudinal trapping frequency uj z « 2tt • 60 Hz. We 



obtain Az 



2X 



AJ Z < 125nm for the position fluc- 



muj^d 

tuations leading to an extra noise of the same order as 
the shot noise limit £^ = OdB for 2000 atoms. Since 
these technical induced fluctuations add to the variance 
caused by the atomic quantum state, their magnitude has 
to be much smaller than the shot noise level in order to 
measure a reasonable amount of number squeezing. As 
a figure of merit, position fluctuations of 60 nm between 
different experimental realizations limit the best observ- 
able number squeezing to £jy ~ — 6dB. This corresponds 
to roughly one percent of the double well spacing in this 
experiment [59 . 



D. Finite temperature. 

The influence of finite temperature is very different in 
the two settings discussed here. While the effective tem- 
perature in the spin degree of freedom is nearly zero in 
the internal case, finite temperature effects easily domi- 
nate the fluctuations when dealing with a BEC in a dou- 
ble well trap |114[ 1122] . The difference emerges from the 
very different initial energy gap between the two modes 
comprising the pseudo spin. When using two hyperfine 
states their energy separation is typically in the order 
of several gigahertz corresponding to temperatures of a 
few hundred millikelvin, much larger than the temper- 
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FIG. 13. Adiabatic versus isothermal spin squeezing in a double well setup. Number squeezing versus mean spin length is 
shown for the same parameters of the experiment reported in 59, 90, 91 (% = 2tt • 0.5 Hz, N = 1600, uj z ~ 2tt • 65 Hz). The 
solid thick black lines are the limits for number squeezing and coherent spin squeezing. The dashed thick black line corresponds 
to the ground state of the Hamiltonian (16). Light gray dashed lines are isothermals at different temperatures indicated on 
the left. The dark gray solid lines correspond to an adiabatic evolution of the system where the number of initially thermal 
populated many-body modes is given in brackets (corresponding to 60 nK, 5nK and 1 nK initially). The regime parameter A 
increases for each line from the right to the left. 



ature in the BEC of a few tens of nanokelvin. Thus, 
using almost noise free coupling fields and short coupling 
pulses much faster than any thermalization timescale, 
spin states with almost zero temperature can be realized 
(see references |111( 1112] for the effect of finite tempera- 
ture in this system.). 



The situation is very different in the external double 
well potential, where the pseudo spin in comprised of the 
two center of mass modes in the two potential wells. In 
a unbiased (5 = 0) situation, there is no or a negligi- 
bly small energy gap separating the two modes such that 
several collective modes are excited even at nanokelvin 
temperatures. This leads to non pure initial spin states 
and to increased spin fluctuations as compared to a coher- 
ent spin state. Nevertheless the adiabatic spin squeezing 
technique results in an effective cooling of the spin fluc- 
tuations such that coherent spin squeezed states can be 
realized [59] . Figure 13 shows the number fluctuations 



and mean spin lengths for different initial temperatures 
comparing the adiabatic method to the isothermal case. 
Parameters are the same as for the experiments reported 
in [59]. Note that the fluctuations for the ground state 
would be orders of magnitude better for the same Hamil- 
tonian parameters, showing the dominating effects of fi- 
nite temperature. 



V. EXAMPLES OF EXPERIMENTALLY 
GENERATED COHERENT SPIN SQUEEZING. 

In the previous sections we took two settings as ex- 
amples for the generation of adiabatic and non-adiabatic 
spin squeezing in a BEC. We will briefly discuss some 
results of these experiments [59, 60 in the following. A 
more detailed description of these experiments can be 
found in EDI EH- 



A. 



Adiabatic spin squeezing in a double well 
potential. 



The experiments discussed in this section have been 
originally published in [59] and in |123j . Approximately 
2000 condensed Rubidium atoms were loaded into a sym- 
metric double trap where initially the potential barrier 
separating the two modes was small. For the initial pa- 
rameters the system was already in the Josephson regime 
due to the rather strong nonlinearity. Next the poten- 
tial barrier was increased adiabatically, driving the sys- 
tem deeper into the Josephson regime. The decrease of 
the tunnel coupling resulted in a decrease of the relative 
number fluctuations directly corresponding to the spin 
fluctuations in J z direction. The fluctuations were redis- 
tributed to the orthogonal spin direction such that the 
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FIG. 14. Number difference fluctuations and mean spin length. The figure shows the result of 1500 measurements with 
experimental parameters chosen to minimize fluctuations in the relative phase (left column) or in the atom number difference 
(right column), a Number difference fluctuations. The solid black line is a Gaussian fit to the left histogram. Its width is 
compatible with the shot noise noise limit taking detection noise into account. For comparison the solid black curve is copied to 
the right histogram whose width is approximately a factor of two smaller, b The lower panel shows a vertical lineup of several 
measurements of the relative phase taken for the same experimental parameters as the number fluctuation measurements above. 
Each row is a horizontal profile of an independent measurement and the gray shading is proportional to the 1-dimensional atomic 
density. Below the ensemble average is shown revealing the mean spin length from the pattern visibility. Small variance in the 
atom number difference coincides with large variance in the relative phase (reduced mean spin length) and vice versa. This 
figure was originally published in 



mean spin length decreased. The simultaneous change 
in the number fluctuations and mean spin length can be 
seen directly on the raw data (figure 14). 



Predicting the optimal time required for an adiabatic 
change of the parameters is not trivial in such an ex- 
periment, especially because additional aspects such as 
particle losses or heating limit the duration of the exper- 
iments. Therefore, the number squeezing was measured 
versus ramp up time of the double well potential while the 
final barrier height was kept constant. Figure [15] shows 
the results of this experiment. In order to explain the de- 
pendence of number squeezing versus ramp up time both 
finite temperature and finite fluctuations of the trap po- 
sition had to be taken into account. They dominate the 
observed fluctuations (details of this experiment can be 



found in [59ll90ll9TirT23]y 

Despite of these limitations coherent spin squeezed 
states in the double well trap were realized with a maxi- 
mal squeezing factor of £| = — 2.3dB [59]. However, the 
implementation of a full atom interferometer beyond the 
standard quantum limit remains an experimental chal- 
lenge with such a setup even though the feasibility of the 
production of coherent spin squeezed states was demon- 
strated. The analog of electromagnetic coupling pulses 
rotating the spin state non-adiabatically needs to be im- 
plemented in order to realize a Ramsey sequence with 
the spin squeezed state as input state. The required non- 
adiabatic changes are difficult to realize without the pos- 
sibility to change the nonlinearity directly (which could 
be done by using a Feshbach resonance [98]) since fast 
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FIG. 15. Timescales for adiabatic barrier ramps. To find the 
optimum ramp times the barrier was ramped linearly from 
a fixed start to a fixed end value with different speeds. The 
optimum is found where number squeezing starts to saturate. 
The hatched area, the theoretical two mode prediction for 
typical temperatures between 20 nK and 40 nK, cannot ex- 
plain the observations. Including position noise (gray area) 
the measured data is reproduced by the theory for fluctuations 

of Az«80nm EgaE]. 



changes in the external trapping potential lead to un- 
wanted excitation of the condensate. 
Using internal hyperfine states to define the two modes 
instead of the center of mass modes these limitations can 
be overcome as presented in the following section. 



B. Non-adiabatic squeezing based on internal 
atomic states. 

When working with internal hyperfine states of the 
atoms fast and precise spin rotations are readily avail- 
able using electromagnetic coupling. The experimen- 
tal challenge is to obtain a sufficiently high nonlinearity 
such that spin squeezed states can be realized within the 
experimental coherence time. Several ideas to increase 
the nonlinear interaction have been experimentally real- 
ized reaching from cavity mediated interactions in ther- 
mal gases [53] [54] to state dependent microwave traps 
on atom chips which allow for the control of the spatial 
overlap of the two spin modes [84]. We focus here on 
another method exploiting an interstate Feshbach reso- 
nance between two hyperfine states of Rubidium around 
9 G [60] R Such a Feshbach resonance may be used to 



3 The two states are the F 
hyperfine states. 



l,rriF = 1 and F = 2,rriF 



FIG. 16. Shown is the measured interspecies scattering length 
aab(B) around the Feshbach resonance [90j IM] . The black 
solid line is a fit to the real part of a a b(B) — a a f(l — 
( B -B 6 -i — 72) ESI revealing an inelastic width of 75 = 4.6 ± 
0.7 mG in accordance with the theoretical prediction [124] . 
The elastic width extracted from the fit is SB — 1.6=b0.2mG. 
The center of the resonance is at Bq = 9.092 G, but the 
absolute magnetic field has an uncertainty of approximately 
10 mG in this measurement [90] [91] • The background scat- 
tering length a a f was chosen as a free parameter and the fit 
reveals a a f = 96.5 db 0.7 as- 



increase the nonlinearity but at the cost of increased par- 
ticle losses [98]. Figure [l6| shows the measured resonantly 
enhanced effective nonlinearity [90l EJ 1107] . The simul- 
taneously enhanced losses limit the available time to gen- 
erate coherently spin squeezed states to a few tens of mil- 
liseconds favoring the fast non-adiabatic "one axis twist- 
ing" method proposed in [10 for the generation of spin 
squeezing [90] [91] II 10] . Coherent spin squeezed states 
with inferred squeezing factors up to £| = — 8.2 dB could 
be realized within an evolution time of 18 ms [60 . The 
use of a spin echo technique and an actively stabilized 
magnetic field [90[ [9TJ I107J was crucial here in order to 
minimize the dephasing effects. 

Given the possibility of precise spin rotations a full 
interferometric sequence can be implemented. In a pro- 
totypical measurement enhanced interferometric sensitiv- 
ity of 15% beyond the standard quantum limit has been 
demonstrated [60] using coherent spin squeezed states 
within the Ramsey interferometer. Figure [IT] shows a 
direct comparison of the Ramsey interferometer perfor- 
mance for coherent spin states as input state versus co- 
herent spin squeezed states. In this direct comparison the 
performance increase due to the spin squeezing is even 
larger (31%). However, it should be made clear that ex- 
perimental noise is included here, that is, the comparison 
is not relative to the standard quantum limit. Measured 
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FIG. 17. Linear versus nonlinear interferometric performance. Shown is the measured atom number difference versus relative 
phase after a Ramsey interferometric sequence. The error bars are two standard deviation statistical errors of the raw data. 
The black (gray) dashed lines are a fit through the upper and lower ends of the two standard deviation error bars for the 
non-linear (linear) interferometer. The horizontal width of the included areas measure the phase precision. The phase precision 
for the non- linear interferometer is increased by 31% as compared to the linear one. 



is the increase in performance due to spin squeezing in 
this particular interferometric setup. 

Several other experiments recently succeeded in the re- 
alization of quantum enhanced interferometers. An in- 
terferometer with 1.1 dB precision gain has been real- 
ized using a cold sample of Cesium atoms where coher- 
ent spin squeezing has been generated by quantum non- 
demolition measurements [85]. Using cavity mediated 
interactions in a cold Rubidium gas to generate coherent 
spin squeezing the time dependent performance of a sim- 
ilar interferometer has been characterized in terms of the 
Allan variance [55] . 

These different experiments show that atomic coherent 
spin squeezing can nowadays be realized with a variety 
of techniques and in quite different systems. It remains 
an open challenge to realize a large amount of coherent 



spin squeezing with atom numbers comparable to those 
used in modern linear interferometers. However, the ba- 
sic concept of spin squeezing and its experimental feasi- 
bility has been demonstrated such spin squeezing might 
help to push the limit of precision measurements in the 
future. 
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